INDEPENDENT PARTICLES IN A DYNAMICAL 
RANDOM ENVIRONMENT 



MATHEW JOSEPH, FIR AS RASSOUL-AGHA, AND TIMO SEPPALAINEN 

Abstract. We study the motion of independent particles in a dynamical random en- 
vironment on the integer lattice. The environment has a product distribution. For the 
multidimensional case, we characterize the class of spatially ergodic invariant measures. 
These invariant distributions are mixtures of inhomogeneous Poisson product measures 
that depend on the past of the environment, and we also investigate the correlations in 
this measure. For dimensions one and two, we also prove convergence to equilibrium from 
spatially ergodic initial distributions. In the one-dimensional situation we study fluctu- 
ations of net current seen by an observer traveling at a deterministic speed. When this 
current is centered by its quenched mean its limit distributions are the same as for classical 
independent particles. 



1. Introduction and results 

This paper studies particles that move on the integer lattice Z d . Particles interact 
through a common environment that specifies their transition probabilities in space and 
time. The environment is picked randomly at the outset and fixed for all time. Given 
the environment, particles evolve independently, governed by the transition probabilities 
specified by the environment. 

We have two types of results. First we characterize those invariant distributions for the 
particle process that satisfy a spatial translation invariance. These turn out to be mixtures 
of inhomogeneous Poisson product measures that depend on the past of the environment. 
Poisson is expected, in view of the classical result that a system of independent random 
walks has a homogeneous Poisson invariant distribution [4, Section VIII. 5]. For d = 1, 2, we 
use coupling ideas from [6] (as presented in [13]) to prove convergence to this equilibrium 
from spatially invariant initial distributions. 

In the one-dimensional case we study fluctuations of particle current seen by an observer 
moving at the characteristic speed. In the present setting the characteristic speed is simply 
the mean speed v of the particles. More generally, the characteristic speed is the derivative 
H'{p) of the flux H as a function of particle density p. The flux H(p) is the mean rate of 
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flow across a fixed bond of the lattice when the system is stationary with density p. For 
independent particles H(p) = vp. 

It is expected, and supported by known rigorous results, that the current fluctuations 
are of order n 1 / 4 with Gaussian limits if the macroscopic flux H is linear, and of order 
n 1 / 3 with Tracy- Widom type limits if the flux H is strictly convex or concave. In statistical 
physical terminology, the former is the Edwards- Wilkinson (EW) universality class, and the 
latter the Kardar-Parisi-Zhang (KPZ) universality class. (See [2] for the physics perspective 
on these matters, and [3, 15] for mathematical reviews.) Our motivation is to investigate 
the effect of a random environment in the EW class. We find that, when the current is 
centered by its quenched mean and the environment is averaged out, the fluctuation picture 
in the dynamical environment is the same as that for classical independent random walks 
[7, 14]. Consistent with EW universality, the current fluctuations have magnitude i 1 / 4 and 
occur on a spatial scale of t 1 ^ 2 where t denotes the macroscopic time variable. However, 
there is an interesting contrast with the case of static environment investigated in [10]. In 
the static environment, the quenched mean of the current has fluctuations of magnitude 
i 1 / 2 and follows a Brownian motion. Preliminary calculations, not yet reported in this 
paper, suggest that under a dynamic environment the quenched mean of the current has 
fluctuations of magnitude t 1//4 and obeys a fractional Brownian motion when the particle 
system is stationary in time. 

We turn to a description of the process and then the results. 

The particles move in a space-time environment uj = (^s,x)( s .x)eZxZ d indexed by a discrete 
time variable s and a discrete space variable x. The environment at space-time point 
(x,s) G 1i d x Z is a vector uj s>x = (oj SiX (z) ■ z G Z d , \z\ < R) of jump probabilities that 
satisfy 

(1.1) < lj 8>x (z) < 1 and ^ u s ,x(z) = 1- 

zGZ d :\z\<R 

R is a fixed finite constant that specifies the range of jumps. From a space-time point (x, s) 
admissible jumps are to points (y, s + 1) such that \y — x\ < R. In the environment u the 
transition probabilities governing the motion of a Z d - valued walk X, = (X s ) s& z + are 

(1-2) P w [X s+1 =y\X s = x]=7r^ s+1 (x,y) = uj SjX (y-x). 

P u is the quenched probability measure on the path space of the walk X, . The environment 
is "dynamical" because at each time s the particle sees a new environment u> s = {^s,x '■ x £ 
Z d ). 

(p,, &) denotes the space of environments u satisfying the above assumptions, endowed 
with the product <r-algebra &. The environment restricted to levels s G {m, ...,n} is 
denoted by 

w m , n = (u s )m<s<n = (w s ,i :m<s<n,x£Z d ). 

Environments at levels generate u-algebras G m ,n = ff{w m ,n}' In these formulations m = 
—oo or n = oo are also possible. T s ^ x is the shift on Q, that is (T StX uj)t t y = u s +t,x+y 
Let P be a probability measure on f2 such that 

(1.3) the probability vectors (w s ,x)( s ,x)eZxZ d are i-i-d- under P. 
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We make two nondegeneracy assumptions. The first one guarantees that the quenched 
walk is not degenerate: 

(1.4) F{3z G 1 d : < u 0fi (z) < 1} > 0. 

Denote the mean transition kernel by piu) = Eir StS+ i(x, x + u). The second key assumption 
is that 

, . there does not exist x G Z d and an additive subgroup 

^ ' ' GCZ d such that £ zeG p(x + z) = l. 

Another way to state assumption (1.5) is that the averaged walk has span 1, or that it is 
aperiodic in Spitzer's [17] terminology. 

To create a system of particles, let {X U ' J : u G Z d , j G N} denote a collection of random 
walks on Z d such that walk starts at site u: Xq' j = it. When the environment u> is 
fixed, we use P u to denote the joint quenched measure of the walks {A A _ MJ }. Under P u 
these walks move independently on Z d and each walk obeys transitions (1.2). 

Further, assume given an initial configuration 7] = (t](u)) ueZ d of occupation variables. 
Variable rj(u) G Z + specifies the number of particles initially at site u. P^ denotes the 

quenched distribution of the walks {X^ ,:> : m £ Z d , 1 < j < t](u)}. Occupation variables for 
all times s G Z + are then defined by 

??(u) 

V'V) =EE ^ x s ,j = x h M £Z + x 7L d . 

When the initial configuration 77 = rjo has probability distribution fi we write Pff = 
JP^(d V ). 

When the environment is averaged over we drop the superscript uj: for any event A 
that involves the walks and occupation variables, and any event B C Q, P^(A x B) = 
J B P^(A)F(ckj). 

It will be convenient to allow the initial distribution [i to depend on the environment: 
fi = fi u . But then it will always be the case that fi u depends only on the past ca-oo-i of the 
environment. Consequently the initial distribution ^ and the quenched distribution of the 
walks P U ({X* '"'} G • ) are independent under the product measure P on the environment. 
And always under a fixed lj, the initial occupation variables {r]o( x )} are independent of the 
walks {Xf > j }. 

The first result describes the invariant distributions of the occupation process ijt = 
(i]t( x )) X £Z d ■ The starting point is an invariant distribution for the environment process 
seen by a tagged particle: this is the process T n ^x n 0J where X, denotes a walk that starts at 
the origin. A familiar martingale argument (Proposition 3.1 in Section 3 below) shows the 
existence of an S^oo.^i-measurable density function / on f2 such that E(/) = 1, K(/ 2 ) < 00, 
and the probability measure Poo(dci;) = f{ui) P(dcj) is invariant for the Markov chain T n) x n 0J- 

For < A < 00 let r A denote the mean A Poisson distribution on Z + . For < p < 00 
and u) G SI define the following inhomogeneous Poisson product probability distribution on 
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particle configurations r/ = (T](x)) x& x d: 

(1.6) ^(dv) = ® T pf ( To >*^(dr](x)). 
Define the averaged measure by 

(1.7) p p = [ n p '"F(du). 



Theorem 1.1. Let the dimension d > 1. Consider independent particles on 7L d in an 
i.i.d. space-time environment (as indicated in (1.3)), with bounded jumps, under assump- 
tions (1.4) and (1.5). 

(a) For each < p < oo, p p is the unique invariant distribution for the process rj. that is 
also invariant and ergodic under spatial translations and has mean occupation j rj(x) dp p = 
p. Furthermore, the tail a-field of the state space ll£ is trivial under p p , and under the 
path measure P^p the process n. is ergodic under time shifts. 

(b) Suppose d = 1 or d = 2. Let is be a probability distribution on that is stationary 
and ergodic under spatial translations and has mean occupation p = J r)(x) dv . Then if 
v is the initial distribution for the process r\., the process converges in distribution to the 
invariant distribution with density p: Pu{vt G •} => P p as t — > 00. 

Under the invariant distribution p p the covariance of the occupation variables is 

Co^fo(0),T/(m)] = p 2 Cov[/H, /(T 0jm w)] = p 2 E[f(oj)f(T , m co)} - p 2 , m e Z d . 

The first equality above comes from the structure of p p : given cj, the occupation variables 
are independent with means E p ' [rj(m)] = pf(To^ m uj). Our next result describes these 
correlations. First some definitions. 

Two auxiliary Markov transitions q and q on Z d play important roles throughout much 
of the paper: 

Q{x,y) = ^2 E [ w o,o (z)u ,x(z + y)} 

(11, 

:E[u 0t o(z)uj 0fi (z + y)\ x = 0, ye 
i p(z)p(z + y - x) 2 ^0, y£ 
and 

(1.9) q(x, y) = q(0, y - x) = ^ p(z)p(z + y-x). 

z&, d 

Think of q as a symmetric random walk whose transition probability is perturbed at the 
origin, and of q as the corresponding unperturbed homogeneous walk. 
For 9 e T d = ( — 7r, 7r] define characteristic functions 

(1.10) f(9) = ^ W0)0 (z)e fc , 

z 

(1.11) \(6)= ^2q(O,z)e ie - z =E\r(0)\ 2 
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and 

(1.12) A(0) = Q(Q,z)e ie ' z = \W(0)?- 

Assumption (1.5) implies that the random walk q is not supported on a subgroup smaller 
than Z d , hence A(0) < 1 for G T d \ {0} [17, p. 67, T7.1]. Define a constant /3 by 

(1.13) = 7^3 f ^—^dO. 
y ' P (2vr) d J Jd 1 - A(0) 

The distribution (7(0,2) is not degenerate by assumption (1.4) and hence A(0) is not iden- 
tically 1. Since also A(0) < \(9), we see that j3 G (0, 1] is well-defined. 

Theorem 1.2. Letd>l. Form€Z d \ {0} 

(1.14) Cov[/ M, f(T , m u)] = - JLL ! cos(0 • m) 1 —^- d9 

(2ir) a Jjd 1 - A(0) 

and 

(1.15) Var[/(w)] = /T 1 - 1. 



The compact analytic formulas (1.13) and (1.14) arise from probabilistic formulas that 
involve the transitions q and q and the potential kernel of q. The probabilistic arguments 
are somewhat different in the recurrent (d < 2) and transient (d > 3) cases. The reader 
can find these in Section 4. 

By computing the integral in (1.14) an interesting special case arises: 

Corollary 1.3. For the simplest case where d = 1 and p(x) + p(x + 1) = 1 for some 
x G Z, the fixed time occupation variables in the stationary process are uncorrelated: 
Cov[/H, f(T , m u)] = /or m / 0. 

We turn to study particle current and for this we restrict to dimension d = 1. Define the 
mean and variance of the averaged walk by 

(1.16) v = xp(x) and <r 2 = x 2 p(x) — v 2 . 

x&L x£Z 

For t G R+ = [0, oo) and r G K, let 

%(^) r?o(aO 

(1.17) r n (t, r) = ^ ^ < [mrtj + ry^} - £ E X W > LnutJ + ry^}. 

x>0 j=l x<0 j=l 

Y n (t,r) represents the net right-to-left current of particles seen by a moving observer who 
starts at the origin and travels to \nvt\ + |_ r \/^J m time [nt\ • 

We look at the current under the following assumptions. Given cj, initial occupation 
variables obey a product measure that may depend on the past of the environment, but so 
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that shifts are respected. Precisely, 

given the environment lu, initial occupation variables (??o( a; ))xez have 
^ 28) distribution ^{drjo) = <g> /j,^(dr]Q(x)) where ^ is allowed to depend 

measurably on u)_oo,-i- Furthermore, n% = Hq°' xU1 ■ 

Let P u denote the quenched distribution P^L of initial occupation variables and walks, and 
P = P w (-)P(du;) the distribution over everything: particles, walks and environments. 
Make this moment assumption: 

(1.19) E[ m (0) 2 ] < oo. 
Parameters that appear in the results are 

(1.20) p =E[ m (x)} and a 2 = E[Var w ( % (0))]. 

Next we describe the limiting process. Let ip v 2{x) = (2ttv 2 )~ 1 ^ 2 exp(— x 2 /2u 2 ) denote the 
centered Gaussian density with variance v 2 , and & u 2(x) = f^ OQ f u 2(y)dy the distribution 
function. Let If be a two-parameter Brownian motion on R + x R and B a two-sided 
one-parameter Browian motion on R. W and B are independent. Define the process Z by 



z i t , r ) = VPo Pa2(t- S )i r ~ x)dW(s,x) 

(121) JJ[o,t]xR 

+ cr / (p a 2 t (r — x)B(x) dx. 
Jr 

{Z(t,r) : t G R + ,r G R} is a mean zero Gaussian process. Its covariance can be expressed 
as follows: with 

(1.22) *„ 2 (a0 = Av(z) - x(l - $„ a (x)) 

define two covariance functions on (R + x R) x (R + x R) by 

(1-23) MM, M) = *^ (t+s) (r - q) - *o*\P-.\(r - q) 
and 

(1.24) r 2 ((s,q),(t,r)) =V a 2 s (-q) + y aH (r)-y a 2 {t+s) (r-q). 
Then 

(1.25) E[Z(s,q)Z(t,r)] = ^((s, g), (t,r)) + a 2 T 2 ((s,q), (t,r)) . 

(Boldface P and E denote generic probabilities and expectations not connected with the 
RWRE model.) 

The theorem we state is for the finite-dimensional distributions of the current process, 
scaled and centered by its quenched mean: 

Y n (t, r) = n~^{Y n (t, r) - [Y n (t, r)]}. 

Fix any N G N, time points < t\ < t 2 < ■ • ■ < tjy G R+, space points r\, r 2 , ■ ■ ■ , rjv £ R 
and an N- vector = (0\, . . . , On) £ WL N . Form the linear combinations 

N N 

Y n (0) = ^O l Y n (t i ,r i ) and Z(0) = £ 0iZ{t u n). 

i=l i=l 
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Theorem 1.4. Consider independent particles on Z in an i.i.d. space-time environment with 
bounded jumps, under assumptions (1.3) and (1.5). Let the lo- dependent initial distribution 
satisfy (1.18) and (1.19). With definitions as above, quenched characteristic functions 
converge in L 1 (F): 

(1.26) lim E|£ w (e iF "W) - B(e iZ ^)\ = 0. 

In particular, under the averaged distribution P, convergence in distribution holds for the 
R N -valued vectors as n — > oo: 

(Y n (t 1 ,r 1 ),Y n (t 2 ,r 2 ), - ■ ■ ,Y n (t N ,r N )) (Z(h, n), Z(t 2 , r 2 ), • • • ,Z(t N ,r N )). 

While we do not have a queched limit (convergence of distributions under a fixed ui), 
limit (1.26) does imply that, if a quenched limit exists, it is the same as we have found. 

Further notational conventions. Denote by N = {1, 2, 3, . . . } and Z + = {0, 1,2,...}. 
Multistep transitions are 

ui,...,u t _ s _iez 

Let S denote the set of all measures fi on (Z + ) zd which are invariant under spatial trans- 
lations. Let S e denote the subset of S consisting of ergodic measures. Let X denote the 
set of measures which are invariant, that is fit = fiS(t) = ft for all t £ Z + (fit, fiS(t) here 
denotes the measure on configurations at time t when the initial measure on configurations 
is fi). E u , E,E V ,~E etc will denote expectations with respect to P, P V ,P etc. Constants 
C can change from line to line. 

2. Coupled Process 

This section describes the coupling that will be used to prove Theorem 1.1. We couple 
two processes rft and Q so that matched particles move together forever, while unmatched 
particles move independently. To do this precisely, choose for each space-time point (x, t) a 
collection Hf )X = { v txi v tx-> v tx : J ^ ^} °f i-i-d. Z d -valued jump vectors from distribution 
u>t x- Given initial configurations rjo and £o> perform the following actions. At each site x 
set 

Co (a;) = Vo(x) A ( (x), ffi(x) = (rj Q (x) ~ Co(x)) + , and /3^(x) = (r/ (x) - ( (x))~ ■ 

£,o(x) is the number of matched particles, while f3 Q (x) count the unmatched + and — 
particles. Move particles from each site x as follows: the £o(x) matched particles jump to 
locations x + Vq' 3 x for j = 1, . . . , £0(^)5 the /3q(x) + particles jump to locations x + Vq£ for 

j = 1, . . . , /3q(x), and the /3q(x) — particles jump to locations x+v £ for j = 1, . . . , Pq{x). 
After all jumps from all sites have been executed, match as many pairs of + and — particles 
at the same site as possible. This means that a H — pair together at the same site merges to 
create a single ^-particle at the same site. Since particles are not labeled, it is immaterial 
which particular + particle merges with a particular — particle. When this is complete we 
have defined the state ($,i(x), f3f(x), /Jf (x)) x&Z d at time t = 1. Then repeat, utilizing the 
jump variables for time t = 1. And so on. 
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This produces a joint process (£t,/3 t , /3 t ) sucn that 

60) = *7t(») A Ct(x), f3f{x) = (r] t (x) - Ct(x)) + , and fc(x) = (r? t (x) - Ctfa)) - - 

The definition has the effect that a matched pair of r\ and £ particles stays forever together, 
while a pair of + and — particles together at a site annihilate each other and turn into 
a matched pair. If we are only interested in the evolution of the discrepancies (/3 t + ,/3 t ~) 
we can discard all matched pairs as soon as they arise, and simply consider independently 
evolving + and — particles that annihilate each other upon meeting. 

If we denote by H = : t G Z, x G Z d } the collection of jump variables, and by Gqj 
the function that constructs the values at the origin at time t: 

then it is clear that the values at other sites x are constructed by applying this same function 
to shifted input: 

(2.1) (Zt{x),0+(x),Pt(x))=G o ,t(6 x Ti o ,6 x (D,O x E). 

Here 9 X is a spatial shift: (O x r])(y) = rj(x + y) and (9 x E) tt y = E t:X + y for x,y G Jj d . In 
particular, if the initial distribution fl of the pair (770? Co) is invariant and ergodic under 
the shifts 9 X , while {Et,x : i G Z,a; G Z d } are i.i.d. and independent of (/7o 5 Co)> it follows 
first that the triple (770, Co>^) is ergodic, and then from (2.1) that for each fixed t the 
configuration (£t, f3^~ , f3f~) is invariant and ergodic under the shifts 6 X . 

Let S, resp. S e , denote the set of spatially invariant resp. ergodic probability distributions 
on pairs (77, £) of configurations of occupation variables. 

Lemma 2.1. Let fl G S. The expectations En\^[x)\ and E^[(3f~(x)] are independent of x 
and nonincreasing in t. 

Proof. The independence of x is the shift-invariance from (2.1). Discrepancy particles are 
not created, only annihilated, hence the nonincreasingness in t. □ 

Proposition 2.2. Let d = 1 or 2. Suppose p G S e . Let £^[r/(0)] = p\ and £^[((0)] = p 2 . 
If Pi > P2, we have 

Efi\ft(0)] = Ep[( m (0) - 0(0))-] as t -> 00. 

Proof. We already know from Lemma 2.1 that E^f3^(0) cannot increase. To get a contra- 
diction let us assume that E^f3f~(0) > 5 for all t where 5 > 0. Since E^/3 t + (0) — Efi(3f~(0) = 
Pi — P2 > 0, we also have Ep,(3^~(0) > 5 for all t. 

At time 0, assign labels separately to the + and — particles from some countable label 
set J and denote the locations of these particles by {u)j~(t), w~(t) : j G J}. Each +/— 

particle retains its label throughout its lifetime. The lifetime of + particle is 

Tj~ = inf{t > : w~j is annihilated by a — particle}. 
If Tj~ = 00, then is immortal. Similarly define rj . Let 

/3±(x) = ^l{ U ; J ± (0) = x,rf > t} 
j 
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denote the number of ± particles initially at site x that live past time t. We would like 
to claim that for a fixed t the configuration {(/^(x), f3^ t (x)) : x £ 7j d } is invariant and 
ergodic under the spatial shifts 9 X . This will be true if the evolution is given by a mapping 
Fqj so that ((3q t {x) , (3$ t (x)) = Fo : t(9 x Vo,@xCo,6xE) for all x G 7h d . Such a mapping can be 
created by specifying precise rules for the movement and annihilation of + and — particles 
that are naturally invariant under shifts. 
Then the ergodic theorem implies that 

E^O) = re lim -±- d £ a.s. 

\x\<n 

Here |x| is the £°° norm: for a vector x = (x\, . . . ,Xd), |x| = maxi<j<^ \xi\. Since particles 
take jumps of magnitude at most R, 

5 < E a Wf(0)] = lim l —— V Bf(x) 

- v n n _^ O0 2n + 1 )d w K J 

\x\<n 

v ' \x\<n+Rt 

The initial occupation numbers of immortal +/— particles are 

The limit exists by monotonicity. This limit produces again a functional relationship of the 
type (2.1): 

(/#~(aO,/W*)) = Km = ^mi^MoACoAH) 

= Fo :OQ (6 x r]o,6 x (,o,6 x E). 

Thereby {{Pq 00 (x),Pq 00 (x)) : x G Z d } is spatially invariant and ergodic. 
By the ergodic theorem again 

^oo(O) = Hm — E ^oo(-) a.s. 

|x|<Tl 

while by the monotone convergence theorem 

^00(0)= lim £ A /3± (0) > 5. 

We have shown that the assumption E^f3^{§) > 5 leads to the existence of positive 
densities of immortal + and — particles. However, a situation like this will never arise 
for d = 1 or 2, the reason being that any two particles on the lattice will meet each other 
infinitely often. More precisely, fix any two particles and let X + and X~ denote the walks 
undertaken by these two particles. Then X + and X~ are two independent walks in a 
common environment oj. Let Yt = — Xf . If we average out the environment, then Yt 
is a Markov chain on Z rf with transition q(x,y) given by (1.8). Our assumptions guarantee 
that this Markov chain is recurrent when d = 1 or 2. Hence Yt = infinitely often as 
claimed. We have thus arrived at a contradiction and the proposition is proved. □ 
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3. Invariant measures 

In this section we prove Theorem 1.1. We begin by deriving the well-known invariant 
density for the environment process seen by a single tagged particle. 

Proposition 3.1. There exists a function < / < oo on SI such that Ef = 1, E(/ 2 ) < oo, 
f(oj) is a function o/£D_oo._i, and 

(3.1) f{oo)= f{T-i,xw)Tr-ifi(x, 0) F- almost surely. 
Proof. For N G Z + define 

(3.2) f N (u)= £>^v |0 M). 

/jv(cj) is S-Ar^i-measurable and a martingale with E/tv = 1- By the martingale conver- 
gence theorem we can define 

f{oj) = lim /at(w) (P-almost sure limit). 

A"->00 

Property (3.1) follows because all the sums involved are finite: 

V/(T_ M ^ 10 (x,0)= lim VMT-x^tt^oM) 

= J im E 7r -^-i--i(^ cc ) 7r -i,o(^ )= E^v-i-oM) 

A — >oa ' — ' ' ' V— >oo ' — ' 

= Jim /v+i(w) = /(w). 

A->oo 

In Lemma 3.3 below we show the L 2 boundedness of the sequence {/at}- This implies that 
/at — > f also in L 2 and thereby implies the remaining statements E/ = 1 and E(/ 2 ) < 
oo. □ 

To prove the L 2 estimate for /jy we develop a Green function bound for the Markov 
chain defined as the difference of two walks. Let Xf and X\ be two independent walks in 
a common environment oj, started at x,y S and = Xf — X\. Under the averaged 
measure Y t is a Markov chain on Z d with transition probabilities q(x, y) defined by (1.8). Y t 
can be thought of as a symmetric random walk on 7L d whose transition has been perturbed at 
the origin. The corresponding homogeneous, unperturbed random walk is Yt with transition 
probability q in (1.9). Write P x and P x for the path probabilities of Y. and Y,. Define hitting 
times of for both walks Yt and Yt by 

(3.3) t = inf{n > 1 : Y n = 0} and f = inf{n > 1 : Y n = 0}. 
Lemma 3.2. There exists a constant C < oo such that for all i£Z d and N e N, 

TV Af 

^g fc (x,0)<C^g fc ( 2; ,0). 

fc=0 fc=0 
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Proof. Suppose we had the bound for x = 0. Then it follows for 

N N N N 

E « fe ( x > °) = E * [ E 1 {n=o>] = ^ [ E E x {^=o} 

fc=0 fc=0 i=0 fc=t 

TV n— i N n—i 

= Y^P x {r = i)Y J Q k (0, 0)<C7£p x (r = i)E 9*(0, 0) 
t=0 fc=0 t=0 fc=0 

= C7^g fc (x,0) 
fc=0 

It remains to prove the result for x = 0. Let do = and 

<7j + i = inf{?i > crj : Y n = and / for some fc E {<Tj + 1, . . . , n — 1}}. 

These are the successive times of arrivals to following excursions away from 0. Let Wj, 
j > 0, be the durations of the sojourns at 0, in other words 

Y n = iff crj < n < Oj + Wj for some j > 0. 

Sojourns are geometric and independent of the past, so on the event {o~j < oo}, 

Let Jn = max{j > : Oj < N} mark the last sojourn at that started by time N. Then 



N J N 

E [ E l{n = 0}] < E [ E Wj\ =E^o[ l{^ < AT} 

i=o 

-£k(l + Jjv). 



fc=0 i=0 j=0 

1 



1-9(0,0)' 

Assumption (1.4) guarantees that q(0,0) < 1. 

It remains to bound £b(l + Jjv) in terms of Y2k=o 9 fe (0j 0). The key is that once the 
Markov chain Y\~ has left the origin, it follows the same transitions as the homogeneous 
walk Yk until the next visit to 0. For z 7^ let 

K Z N = inf{fe > 1 : T-f + T| + • • • + T% > N} 

where the {Tf} are i.i.d. with common distribution P z {f & • }• Imagine constructing the 
path Yfc so that every step away from is followed by an excursion of Y& that ends at 
(or continues forever if is never reached). The step bound (1.1) implies that Po{|Yi| < 
2R} = 1. Then there is stochastic dominance that gives 

(3.4) E (J N )< Yl E ( R n)- 

z^0:\z\<2R 

By T32.1 in Spitzer [17, p. 378], for z / 

(3.5) hm g— — = a(z) 
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where the potential kernel a is 

n n 

(3.6) a(z) = lim { ]T q k (0, 0) - ]T q k (z, 0)). 

k=0 k=0 

By P30.2 in [17, p. 361] a(z) > for all z ^ for d = 1,2. For d > 3 by transience 

oo oo oo 

a(z) = Y^q k (0,0) -J^M) = (1 - F(z, 0)) £ g fc (0, 0) > 

A;=0 fc=0 fc=0 

where F(z, 0) = P^j Y n = for some n > 1} < 1. 

From (3.5) and a(z) > 0, there exist < c(z), C(z) < oo such that for all n, 

c(z)P (f > n) < P(T? >n)< C{z)P (f > n) 

and hence 

c{z)E [fAN] <E [T{ A N] < C(z)E [f AN]. 
By Wald's identity and some simple bounds (see Exercice 4.1 in [5, Sect. 3.4]) 

£7(Tf A JV) ~ L J ~ S(Tf A N) 
Let {fj} be i.i.d. copies of f from (3.3) and put 

M N = inf{& > 1 : fi + f 2 H h f fc > TV}. 

Then we have a similar relation: 

N = r n 2iV 

< [Mjv] < 



Eq(tAN) - L J ~P (TAAf) 
Combining the above lines: 

(3J) ^ £ £ c(z)£~o [f A N] £ 

Considering excursions of the y-walk away from 0, 

N N 

(3.8) E [M N ] <l + ^ fc (0,0)< 2^(0,0). 

fc=0 fc=0 

The proof is now complete with a combination of (3.4), (3.7) and (3.8). □ 

From the previous lemma follows the 1? estimate for /jy which completes the proof of 
Proposition 3.1. 

Lemma 3.3. There exists a constant C < oo such that E(/jy-) < C for all N. 
Proof. By translations 

Wn) = ^E7r_jv,o(z, 0)7r_tf,o(z,0) 



J>P-{^ = ^ = ]T^( y ,0). 
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By the submartingale property E(/jy) is nondecreasing in N. Hence it suffices to show the 
existence of a constant C such that 

N 

(3.9) ^E(/|)<CiV for all N. 

From above, by Lemma 3.2 and the spatial homogeneity of the Y-walk, 

N N N 

= £ J> fe (*,0) < C7£ 

fc=0 fc=0 x k=0 x 

N 

fc=0 x 

Property (3.1) implies that the probability measure W ) O0 (duj^) — f(uS} ]P((ic^) is invariant 
for the process T n ^x n ^- Recall from (1.6) the product measure 

(3.10) n p,u (drj) = (g) V f ^' u \dri(x)) 

x&Z d 

where T x is Poisson(A) distribution. By the definition of /, fi p,UJ depends on oj only through 
the levels Q-oo-i. 

Lemma 3.4. The following holds for P-a.e. oj. Let tjq be fj, p ' u '-distributed. Then for all 
times t G Z +; under the evolution in the environment oj, rjt is [i p,Tt ' 0UJ -distributed, and in 
particular independent of the environment ojt at level t. 

Proof. Consider the evolution under a fixed u). The claim made in the lemma is true at time 
t = by the construction. Suppose it is true up to time t — 1. Then over x£Z the variables 
r)t-i(x) are independent Poisson variables with means pf(T t -\ jX oj). Each particle at site x 
chooses its next position y independently with probabilities 7r^ 1 t (x, y). As with marking a 
Poisson process with independent coin flips, the consequence is that the numbers of particles 
going from xtoy are independent Poisson variables with means pf(Tt-i tX u>)ir£_i t (x, y), over 
all pairs (x,y). Since sums of independent Poissons are Poisson, the variables (r]t{y))yeZ 
are again independent Poissons and r)t(y) has mean 

X z 

= Pf( T t,yU). 

The last equality is from (3.1). 

We have shown that r]t = ( r lt{y))y£Z d nas distribution /j 1 p,Tt '° UJ . This measure is a function 
of w -oo,t-i> hence under P independent of tot- D 

Recall from (1.7) the averaged measure \x p = J n p,U} F(doj). 

Lemma 3.5. The measure fi p is invariant and ergodic under spatial shifts 6 X . The tail 
a -field of the state space is trivial under fi p . 
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Proof. Invariance under X comes from fx p ' u o 0~ Y = ^/Vo,<b<«> anc j invariance of IP. Er- 
godicity will follow from tail triviality. 

Let B C Z+ d be a tail event. Then by Kolmogorov's 0-1 law ^^(B) E {0, 1} for each oj. 
We need to show that ji p)UJ {B) is P-a.s. constant. For this it suffices to show that p p,UJ {B) 
is itself (almost surely) a tail measurable function of us. 

Consider a ball A = {(s, z) : \s\ + \z\ < M} in the space-time lattice Z rf+1 . Since the step 
size of the walks is bounded by R, for each N > 1 



is a function of the environments {uj S)Z : s < —1, \z — x\ < R\s\}. Consequently, if \x\ > 
(R + 1)M, the entire sequence {/jvC?b x u)} is a function of the environments outside A, 
and then so is (almost surely) the limit f(To x ui). Since B is tail measurable, ^'^(B) is a 
function of {/(Tq^uj) : \x\ > (R + 1)M} and thereby a function of environments outside 
A. Since A was arbitrary, we conclude that fi p:UJ (B) is (almost surely) a tail measurable 
function of u. □ 

Proof of the first part of Theorem 1.1 ( except for uniqueness) . Invariance of \i p for the pro- 
cess follows by averaging out uj in the result of Lemma 3.4. Spatial invariance, ergodicity 
and tail triviality of fi p are in Lemma 3.5. That f ?](0) dfi p = p follows from the definition 
of n p . 

We prove the ergodicity of the process 77. under the time-shift-invariant path measure 
PaP- We use the notation [i p also for the joint measure p p (dui, drf) = ¥(dco)fi p ' ul (dr]) and not 
only for the marginal on r/. Let J be the u-algebra of invariant sets on the state space of 
the particle system: 



By Corollary 5 on p. 97 of [12] it suffices to show that J is trivial. We establish triviality of 
J by showing that E pP [ip | J\ is almost surely a constant for an arbitrary bounded cylinder 
function tp on . 

Let r) a ' b denote the configuration obtained by moving one particle from site a to site b, 
if possible: rj a,b = n if rj(a) = 0, while if 77(a) > 0, 



Lemma 3.6. There exists a version of E pP [ip \ J\ such that for all rj G Z^ and a, b £ Z rf , 



Proof. By Corollary 2 on p. 93 of [12], we can define a version ip of E pP [ijj \ J] pointwise by 




z&Z d 



J = {B C Zf : l B (v) = P v im e B} for u^-a.s. rj } 




E pP [i;\J}{r ] ) = E pP [ij\J}{r ] a > b ). 




n-l 



We show that VK 7 ?) = 4>{f] a ' h )- 
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Assume 77(a) > 0. Consider the basic coupling P„„a,b of two processes (rjt,Ct) with initial 
configurations (770, Co) = ( 7 ],'n a ' b ), as described in Section 2. Let 

a = inf{rj : ip(7] s ) = tp(Cs) for all s > t}. 

We observe that P ri ^a,b{a < 00} = 1 in all dimensions. In dimensions d £ {1,2} the 
irreducible g-random walk is recurrent, hence the two discrepancies of opposite sign that 
start at a and 6 annihilate with probability 1. In dimensions d > 3 the discrepancies 
are marginally genuinely d-dimensional random walks by assumption (1.5). Thus they are 
transient, and so either the discrepancies annihilate or eventually they never return to the 
finite set of sites that determines ip. 

The conclusion of the lemma follows: 

\E r ,[%b{r] t )] - E va , b [^(i ]t )]\ < 2\\ijj\\ 00 P ri>v a, b {a >f}^0 ast^ra. □ 

Lemma 3.7. Suppose h is a bounded measurable function on such that for all a, b G Z d , 
h(rj a,b ) = h(rj) [i p -a.s. Then there exists a tail measurable function hi such that h = h\ 
jjL p -a.s. 

Proof. To show approximate tail measurability we approximate by a cylinder function and 
then move particles far enough one by one. (This trick we learned from [16].) Let rf denote 
the configuration obtained by removing one particle from site a if possible: 

77 «(x) = /^ (a) ~ 1)+ X = a 
1 n(x) x / a. 

Let e > 0. Pick a bounded cylinder function h such that E pP \h — h\ 2 < e 2 . For each 
to pick b(ui) E Z d so that f (Tqu^oj) > 1/4 and h does not depend on the coordinate 
77(6(0;)). Such b(cj) exists a.s. by the ergodic theorem since E/ = 1. Choose b(u) so that it 
is a measurable function. Since h(w) = h{rj a,h ^) fi p (duj,drj)-a.s. and h(rj a ) = h(rj a ' b ^) by 
choice of b(ui), 

f \h( V )-h(r, a )\^(d V )< f l { , (a)>0} |n(77 a ^))- h(r, a ^)\^(du;,dr,) 

+ J i Ma)>0} \Hv a )-Hv a )\^(.dv). 

In the next calculation we bound the first integral after the inequality. Write 77 = 
(77', 77(a), 77(6(0;))) to make the coordinates at a and 6(0;) explicit. Change summation indices 
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and apply Schwarz. 

1{„(«)>0} - h(r) a > b ^)\ » p (dw, dri) 

E ^2 T Pf(T , a u ) )^ T pf(T 0Maj)U )^ J \h(rf,k- 



k>0 
£>0 



h(r)',k — 1,£ + 1)\ (Ji p ' w (drf) 



E \ - f(T 0t gU}) _ U _ T pf(T , a uj) ( )T pf(T 0Mu>) uj) ( x 

^ m+1 /(T , feHW ) ^ lJ 

n>0 



J Xl(y » ! 77(a) + 1 /(T 0j6 ( w )w) 

n>0 

x {E» P \h-~h\ 2 } 1 ' 2 

< Ce. 

The finite constant C above comes from the property f {Tq^^lo) > 1/4 and C is indepen- 
dent of p, h and h. An analogous argument (but easier since we do not need the b(u>)) 
gives 

/ l Ma )>o}\kv a )-h(v a )\» p (dv)<Ce. 

Since e > was arbitrary we have h(rj) = h(rj a ) fj, p -&.s. 

Applying the mapping r] \— > rf repeatedly shows that, for any given finite A C Z d , h 
equals a.s. a function g\ that does not depend on (rj(x) : x G A). As A increases along 
cubes, the limit hi of the g^s is tail measurable. □ 

We can now conclude the proof of ergodicity. Lemmas 3.6 and 3.7 show that E pP [ip \ J] 
is n p -a,.s. tail measurable, and hence a constant by Lemma 3.5. □ 

3.1. Proof of uniqueness. In this subsection, we complete the proof of part (a) of Theo- 
rem 1.1 by showing that fi p is the unique invariant distribution with the stated properties. 
We also prove the second part of the Theorem 1.1. The proof of uniqueness uses standard 
tricks employed in interacting particle systems [8] . We will arrive at the proof of uniqueness 
through a sequence of lemmas. 

For two configurations rj, £ of occupation variables, we say that rj < £ if rj(x) < £(x) for 
all x. For two probability distributions fj,, v on the configuration space, we say /i. < v if 
there exists a probability measure fl on pairs (77, £) of configurations of occupation variables 
such that p(r] < C) = 1, and the marginals are \i and v. For a convex set A, A e will denote 
the set of extremal elements. 
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Recall that S, resp. <S e , denotes the set of spatially invariant resp. ergodic probability 
distributions on pairs (rj, Q of configurations of occupation variables. Let I denote the set 
of invariant probability distributions on pairs of configurations of occupation variables. 

Lemma 3.8. If pi < p2 then p Pl < p p2 . 

Proof. Define a measure p on the product configuration space as 

£{foC)e-}=E[M w {foC)e-}]. 

Here under p u , the occupation variables rj(x) are independent Poisson with means pif(To x w) 
and = n(x) + r y(x), where 7(2:) are independent Poisson with means (p2 — pi)f(To tX u). 
It is easy to see that p has the required marginals. □ 

We state the next two lemmas without proof. The proofs can be found in Lemmas 4.2 - 
4.5 of [1]. 

Lemma 3.9. We have 

(a) If pi, P2 £lfl 5, there is a p G I OS with marginals p\ and p2- 

(b) If p\,p2 G (Xn5) e; there is a p G (I fl 5) e with marginals p\ and p2- 

Lemma 3.10. If p, E (X n S) e and p{(r], () : rj > ( or ( > n} = 1 then 
p{(v, : V > C} = 1 or p{(r], C) : ( > rj} = 1. 
A crucial lemma needed in the proof of uniqueness is the following. 
Lemma 3.11. Let p e S e such that J[r/(0) + C(^)]dp < 00. Fix x / y G Z d . Then 
lim /2t{(7/,C) : r](x) > ((x) and rj(y) < ((y)\ = 

Proof. Our proof employs some of the notation developed in Section 2. Fix a positive integer 
m. Let / = [—m,m] d and let B be the event that I contains both + and — particles. The 
theorem will be proved if we can show that pt{B) — > as t — > 00. So let us assume to the 
contrary that we can find a sequence tk t 00 such that 

(3.11) p tk (B)>5>0. 

By our assumptions on the environment, we can find a positive integer T = T(m) and a 
positive real number p = p(m) > such that 

min P I X x and X v meet by time T ) > p. 

x,y£l V / 

Let A(t, y) denote the event that a + or a — particle present in the cube y + / at time t 
has been annihilated by time t + T. It is clear that 

(3.12) P[A{t,y)\fkM >P^B{0 y (vtXt)} a.s. 

For what follows, assume that t n+1 — t n >T. Also let 4>t(x) = /3 t + (x) + /3^~(x) be the total 
number of discrepancy particles at x at time t. Let n = l(2m + 1) + m for a positive integer 
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/ and divide the cube [— n, n] d into cubes of side length 2m + 1. We have 

ye[-n-RT,n+RT} d 
(2l+l) d 

(2n + 1)~ 

where C(j) is the center of cube j. Taking expectations and letting n — > oo, we get by the 
ergodic theorem, 

(2l+l) d 

1 •M(t fcl O(j)) 



1 (2t+lj» 

— M^Cy) ■lA(t lbl C(j)) 



_ (2Z+l) a 

Efi<lH h +T(0) < EfiMO) + Km sup £ J ^ MWfo> CtJ} • I 

3=0 



It follows from (3.12) and (3.11) that 



l B{o C {j){m k Xt k )} ■ ^A(t k ,cu)) > ph k {B) > pS. 



We thus have 

^ +r (o)<^ fc (o)- (2m + 1)d - 

We can conclude from Lemma 2.1 that En&uify ~ > —oo. But this is a contradiction since 
En(f>t(0) > 0. The proof of the lemma is complete. □ 

Lemma 3.12. If p,±, p2 G (In <S) e and £ ft r)(0) < oo /or i = 1, 2, i/ien /ii < p,2 or p,2 < p\. 

Proof. Prom Lemma 3.9, we can find p, G (I n 5) e with marginals p,\ and p,2- Using the 
ergodic decomposition of stationary measures [18, Theorem 6.6], 

P{(V,C) ■ V(z) > C(x) and rj(y) < ((y)} = / v{(rj,C) ■ v(z) > C( x ) and V(v) < CG/)}r(dz>), 

JSc 

for a probability measure r on 5 e . On applying the operator S(t) to both sides of the above 
equation, we observe that the right hand side goes to 0. We thus get 

/HfoC) :r/<CorC<??} = l 

An application of Lemma 3.10 completes the proof. □ 

Proposition 3.13. I/jtiG (In5) e and p = E^O) < oo then p = p p0 . 

Proof. Since p p G 5 e nl, it follows that p p G (Xn«S) e . We can then conclude from Lemma 
3.8 and Lemma 3.12 that there exists a po such that li < p p for p > po and p > p p for p < p$. 
Now fix pi < po < p2- For all (x±, X2, • • • , x n ) G (Z d ) n and all (k\, fe, • • • , k n ) G (Z + ) n , we 
have 

p pl (r](xi) > ki,l < i < n) < p(rj(xi) > fej, 1 < i < n) < p P2 (rj(xi) > < i < n) 

The first inequality (resp. the second inequality) above can be seen by looking at the 
coupled measure p corresponding to p pl (resp. p P2 ) and p so that p{r\ < £) = 1 (resp. 
AC 7 / > C) = !)• Now let p\ t and y»2 4 to see that p has the same finite dimensional 
distributions as p p0 . □ 
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Proof of the remaining parts of Theorem 1.1. We first prove that pP is the unique measure 
with the stated properties in part (a) of Theorem 1.1. Indeed, let p be another measure with 
those properties. Since /i 6 iS e nI, we can conclude that p E (In S) e . From Proposition 
3.13, we must have that p = p p . 

We now turn to part (b) of the theorem. Let v be a probability measure on that 
is stationary and ergodic under spatial translations and has mean occupation j £(0)dv = 
p. Denote the occupation process with initial distribution v by Q- Utilizing the ergodic 
decomposition theorem [18, Theorem 6.6], find p E S e with marginals p p and v. Let pt be 
the time t distribution of the joint process (r]t, Ct) coupled as described in Section 2. 

Initial shift invariance implies that mean occupations are constant p throughout time 
and space: 

E u [Ct(x)}= [E{Y,c(y)<Ay^)}^ d O = T, E Kt(y^)) U{y)dv = P . 

y y 
Chebyshev's inequality and Tychonov's theorem can be used to show that the sequence 
{p t }t&+ is tight. 

Let v be any limit point as t — > oo. Then by Proposition 2.2 C) : V = C} = 1- This 
proves that P u {Ct E •} =^ p> p . This completes the proof of Theorem 1.1. □ 

4. COVARIANCES OF THE INVARIANT MEASURES 

Define the Green's functions for both q and q walks by 

N N 
G N (x,y) = ^2q k (x,y) and G N (x, y) = ^ q k (x, y). 
k=0 k=0 

Additionally, for the q walk define the potential kernel by 

(4.1) a(x)= lim (G N (0,0) -G N (x,0)). 

In the transient case d > 3 the limit above exists trivially, since 

oo oo 

G{x, y) = ^2 q h ( x i y) < °° and G(x, y) = ^2 q k ( x i v) < °°- 

k=0 k=0 

So for d > 3 

(4.2) a(x) = G(0, 0) - G(x, 0). 

For the existence of the limit (4.1) in the recurrent case d E {1, 2} see Tl on p. 352 of [17]. 
In all cases the kernel a(x) satisfies these equations: 

q(0, z)a(z) = 1 and q(x, z)d(z) = a(x) for 

z z 

The constant (3 defined by (1.13) has the alternate representation 

(4.3) = J2q(O,z)a(z). 

z 

We omit the argument for the equality of the two representations of j3. It is a simple version 
of the one given at the end of this section for (1.14). 
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To prove Theorem 1.2 we first verify this proposition and then derive the Fourier repre- 
sentation (1.14). 

Proposition 4.1. Let d > l. For m e Z d \ {0} 

(4.4) Cov[/(w), f(T , m oj)] = /T 1 ?(0, z)[d(-m) - a(z - m)] 

z 

and 

(4.5) Var[/(w)] = /T 1 - 1. 

A few more notations. Recall that Y n denotes the Markov chain with transition q and 
Y n the q random walk. Successive returns to the origin are marked as follows: 

(4.6) r = and for j > 0, tj = inf{n > 7j_i : Y n = 0}. 

Abbreviate r = t\. The corresponding stopping time for Y n is f . For m S Z d and A > 1 
abbreviate 

C N (m) = Cav[fN(oj),fN(T 0>m (jj)] = ^ Cov[7r_7v >0 (z, 0), 7r_Ar i0 (u;, m)] . 

z,w£l d 

Define also the function 

h(y) = Cov[7r ,i(0, y + z), 7r 0j i(0, z)] = g(0, y) - g(0, y), y £ Z d . 

Symmetry h(—y) = h(y) holds. 
Lemma 4.2. In a// dimensions d > 1, 

(4.7) CW(m) = %)Gjv-i(y,m). 

Proof. The case A = 1 follows from a shift of space and time. To do induction on A 
use the Markov property and the additivity of covariance. Abbreviate temporarily K XtV = 
Tt-N-N+iixiV) and recall that the mean kernel is p y - x = Ek X]!/ . 

Cjv(m) = Y Cov[K ZiZl 7r_Ar +li0 (zi,0) , K u , iU , 1 7r_ A r + i i o(wi,m)] 

Z,Zl,W,Wl 

(4.8) = Y {Cov[(k ZiZ1 - p Zl - z )7r^N+ifi(zi , 0) , (Kui,uii - Piui-wK-JV+l.oCltfl, m)] 

Z,Zl,W,Wi 

(4.9) +Cov[jp 2l _«7r_Ar+i,o(^i,0), - p TOl _ w )7r_jv+i,o(t«i, m)] 

(4.10) +Cov[pg 1 _ a 7r_jv+i,o(^i,0) , p Wl -wn-N+i,o(m,™>)] }. 

Working from the bottom up, the terms on line (4.10) add up to C^-i{m). The terms on 
line (4.9) vanish because k W;W1 — p Wl - w is mean zero and independent of the other random 
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variables inside the covariance. On line (4.8) the covariance vanishes unless z = w. Thus 
by rearranging line (4.8) we get 

C N {m) - CW_i(m) = line (4.8) 

= ^2 Cov(K 2j2l ,K 2i „, 1 )E[7r_ A r + i i o(zi,0)7r_ A r + i i0 (^i,m)] 



2,21, Wl 



^ Cov(k ,x, K0,x+y) ^2 E [ n -N+l,o(y^ m + ^-N+lfliO, tj\ 



In the recurrent case we will use Abel summation, hence the next lemma. 
Lemma 4.3. Let d G {1,2}. For x,m G Z d , the limit 



(4.11) 



a(x, m) = lim s k [q k (0, m) — q k (x, mj) 
s ' 1 k=0 



exists. For m = the limit is 
(4.12) 

and for m/0 



a(x, 0) 



a(x) 



(4.13) 



a(x,mj 



a(x) 



q(0, z) [a(— m) — a(z — m)] — a(— m) + a(x — m). 



Proof. Let s vary in (0, 1) and let 



r-l 



T-l 



C/(x,m, S ) = ^[^ S fc l{y fc = m}] -+^[^l{n = m}] = l/(x,m). 
fc=o s fc=0 

Decompose the summation across intervals [tj, Tj+i) and use the Markov property: 

OO 71— 1 OO r j + l — 1 

«V m) = E x [ £ s fc l{y fc = m}] +J2 E * E sfc_Tj X { y * = m > 



fc=0 



fc=0 



j'=i 



k=T; 



U(x, m, s) + ^ J B j; (s T )Eo(s r ) j_1 ^(0, m, s) 

3=1 

I7(x, m, s) + ^^. [/(O, m, s). 

i - -E;o(s r ) 



From this, 



(4.14) ^/(/(O,™)-^,™))- 



k=0 



U(0,m,s) — U(x,m,s). 



We analyze the quantities on the right in (4.14). 



□ 
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Suppose first i/O. Then U(x, m) is the same for the Markov chain Yk as for the random 
walk Yfc because these processes agree until the first visit to 0. In the notation of Spitzer 
[17], with a bar added to refer to the random walk Y/., gio}( x i m ) = U(x,m). By P29.4 on 
p. 355 of [17] and DILI on p. 115 that connects a and A, for recurrent random walk 

U(x,m) = gr \(x,m) = A(x,0) + A(0,m) — A(x,m) = a(x) + a(—m) — a(x — m). 
For x = we have U(0, 0) = 1, and for ra/0, 

JJ(0, m) = ^ y)U (y, m) = ?(°> V) [^(v) + a{-m) - a{y - m)] 

y^O y^O 

= P + XI q ( ' y "> ["(~ m ) - a(y - m)] . 
y 



For the asymptotics of the fraction on the right in (4.14) we can assume again x ^ for 
otherwise the value is 1. It will be convenient to look at the reciprocal. A computation 
gives 

1-£q(Q _ EZoS k Po(r>k) 



Again we can take advantage of known random walk limits because both x, z 7^ so the 
probabilities are the same as those for Yj.. By P32.2 on p. 379 of [17], as s /• 1, for recurrent 
random walk the above converges to (note that E x (r) = 00) 

E, a(z) 8 

z ^ a ( x ) a ( x ) 

Letting s / 1 in (4.14) gives (4.12) and (4.13). □ 

For m = we can obtain the convergence as in (4.1) without the Abel summation. But 
we do not need this for further development. 

Proof of Proposition 4-1- Since /jv — > f in L 2 (F), the covariance in (4.4) is given by the 
limit of Cjv(m), so by (4.7) 

Cov[/(w),/(T 0m a;)] = Km { V q(0, y)G N - 1 (y, m) - V q(0, y)G N - 1 (y, m) \. 

y y 

Next, 

'^2q(0,y)GN~i{y,m) = GN(0,m) - 8 , m = q N (0,m) - 5 ,m + Gjv-i(0,m). 
y 

Since the Markov chain q follows the random walk q away from it is null recurrent for 
d = 1,2 and transient for d > 3. So q N (0,m) — > [9, Theorem 1.8.5]. Thus the limiting 
covariance now has the form 

(4.15) -S , m + hm yVo.j/MGW-iCO.m) -G^_i(y,m)]. 

N~>oo 1 — ' 
2/ 
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At this point the treatment separates into recurrent and transient cases. 
Case 1. de {1,2} 

Convergence in (4.15) implies Abel convergence (Theorem 12.41 in [19] or Theorem 1.33 
in Chapter III of [20]), so the limiting covariance equals 

oo 

-S 0>m + lim V q(0, y) V s k (q k (0, m) - q k {y, m)) . 
y k=0 

By substituting in (4.12) and (4.13) we obtain (4.5) and (4.4). 
Case 2. d > 3 

In the transient case we can pass directly to the limit in (4.15) and obtain 

(4.16) Coy[f(oj),f(T , m uj)} = -6o tm + ^q(0,y)[G(0,m)-G(y,m)]. 

y 

The sum above can be restricted to y ^ 0. By restarting after the first return to 0, 



r-l 



(4.17) 
Next, 



G(y,m) =E y [j2 ^ = m i] + P v( T < °o)G(0,m). 



k=0 



(4.18) 



G(0,m) = ^ J E [l{r i <oo} E l{Y k = m} 

j=0 k=T 3 

OO T — 1 

= E P ^ T < °°) JjE; o [ E 1 { y * = m > 



j=0 



fc=0 



1 



r-l 



Po(r = oo) 



>0,m 



+ (1" )^g(0^)E z [^l{y fc = m} 



2^0 



k=0 



Now consider first m^O. Combining the above, 

Cov[/H, f(T , m cj)\ = E ?(0, y){G(0, m) - G(y, m)} 



E *(°' » ) { p £ = ^ E «(°. ^ [ E w = ™>] -ME ^ = m >] } 



. o(r = 00) 
using equality of q and away from 



E«(°'tf) 

Po(r = 00) 



P y (r = 00) 
Po(r = 00) 



fc=0 



r-l 



k=0 



r-l 



E 9(0, z )e z [ e = "»}] - E y [ E i{n = m>] } 



2^0 
r-i 



fc=0 



fc=0 

f-1 



, 1-^2 e «(°» ^ [ E ^ = m >] - E f [ E = m > 

°^ ' 2^0 fe=0 yj^O k=0 
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applying (4.17) and (4.18) to the q walk 
Po(f = oo) 



Po(t 


= oo) 


Po(f 


= oo) 


Pq(t 


= oo) 


Po(f 


= oo) 


Po(r 


= oo) 



Y,V&z){G(z,m) - P z (f < oo)G(0,m)} - P (f = oc)G(0, 

z^O 

z)t.G(z,m) - P z {f < oo)G(0,m)} 



2^0 



p { L = ™| E ( r = *o 



^(0,z)[<5(z,m)-G(0,m)]. 
To finish this case, note that 

p = E 9(0, *)a(*) = E g(o, *)(G(o, o) - G(z, o)) = Y, «(o, *) pf-I°°i 

Pq(t = oo) 
^o(r = oo)' 

We have arrived at 

Cov[/(w), /(T , m w)] = ^ g(0, z) [a(-m) - a(z - m)} . 

z^0 

Return to (4.16)-(4.18) to cover the case m = 0: 

Cov[/(o;),/H] =J2my)[G(0,0) - G(y,0)} - 1 = ^ q(0,y) ^ = °°\ - 1 

= ^r = oo)_ i = rl _ L 
P (t = oo) 

This completes the proof of Proposition 4.1. □ 

Completion of the proof of Theorem 1.2. It remains to prove the Fourier representation 
(1.14) from (4.4). In several stages symmetry of a and the transitions is used. 

Cov[/(w), f(T , m oj)] = /T 1 ^(0, z)[a(m) - o(m - z)] 

z 

AT 



fc=0 a 
-1 



J im T^EE^M / [e-^^)-e-«^]A*(e)c» 

(27T)' 1 ^ ,/ T d 



fc=0 ^ 
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(2t) 



The last equality comes from < A(0) < 1 for 6 £ T d \ {0} and dominated convergence. 
The ratio (1 — A(#))/(l — \(9)) stays bounded as 9 — > because both transitions q and q 
have zero mean and q has a nonsingular covariance matrix [17, P7 p. 74]. □ 



5. Convergence of centered current fluctuations 

We prove Theorem 1.4 by proving the following proposition. Recall the definition of the 
current Y n (t,r) from (1.17), and let {Z(t,r) : (t,r) G M + x M} be the mean zero Gaussian 
process defined by (1.21) or equivalently through the covariance (1.25). Recall also the 
definitions 

Y n (t, r) = n-^^Ynit, r) - E"[Y n (t, r)]}, 

N N 

Y n (e) = Y,^Yn(ti,n) and Z(G) = £ n). 

i=l i=l 

Proposition 5.1. 

(5.1) E"[exp{iY n (6)}) ->■ E[exp{iZ(0)}] in ¥ -probability. 



The remainder of the section proves this proposition and thereby Theorem 1.4. We write 
Y n (0) as a sum of independent mean zero random variables (under P u ) so that we can 
apply Lindeberg-Feller: 



and 



(5.4) 



v 



i=l 



(5.2) Y n (0) = n- 1 / 4 Y J e i{ Y n(tun)-E u >Y n (t l ,r l )} =W n = £ U„ 
with 

(5.3) U m = ^2d i (u m {t i ,r i )l{m>0}-V m (t l ,r l )l{m<0}), 



N 



i=l 



Vo(m) 

U m (t,r) = n" 1 / 4 £ < Lnt-tJ + r^} 



n" 1 / 4 £ w (r|o(m))P u (IS < [m>*J +r 



% (m) 

F m (i,r) = n" 1 / 4 ]T > [nutj +ry^} 

_ n -i/4^ (r?o(m))p c (x m > LnwtJ +rv^). 
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The variables {U m } m £z are independent under P^ because initial occupation variables and 
walks are independent. We will also use repeatedly this formula, a consequence of the 
independence of tjq and the walks under P^: 



(5.5) 



E w [U m (t,r) 2 ] =n- 1 / 2 Var w l ^ < [nvt\ +r^}) 



= n- l l 2 E u { m {m))P UJ {X™ < [nvt\ +r^i)P UJ {X™ > [nvt\ +r^i) 

+ n- 1/2 Var w (r ?0 (m))P CJ (X^ < [nvt\ + ry/n) 2 

and the corresponding formula for V m (t,r). 

Let a{n) / oo be a sequence that will be determined precisely in the proof. Define the 
finite sum 



(5.6) 



IT", 



Y] u m . 

\m\<a(n)y/n 



We observe that the terms |m| > a(n)y / n can be discarded from (5.2). 
Lemma 5.2. E\W n — W*\ 2 — » as n — > oo. 

Proof. By the mutual independence of occupation variables and walks under P^, and as 
eventually a(n) > \ri\, the task boils down to showing that sums of this type vanish: 

2 ■ 

' m>a{n)y/n m>a(n)y/n 

< ri'^E [ EU (vo(rn))+Vai u (r }0 (m))]P u '{X^ t ' j < [nvtj +r^} 

m>a(n)y/n 

r^Jn — m} 

m>a(n)^n 

' X nt - [nvt\ 



\ 2 1 

£ U m (t,r)\ =E Yl E^[U m (t,r) 2 ] 



CE 



n 



r + a{n) 



Under the averaged measure P the walk X s is a sum of bounded i.i.d. random variables, 
hence by uniform integr ability the last line vanishes as a(n) oo. There is also a term for 
m < a(n)^/n involving V m (t,r) that is handled in the same way. □ 

The limit 6 ■ Z in our goal (5.1) has variance 

(5.7) o- 2 e = £ d^por^^n),^,^)) +a 2 T 2 ({t l ,r i ),(t J ,r J )) 

and the two T-terms, defined earlier in (1.23) and (1.24), have the following expressions in 
terms of a standard 1-dimensional Brownian motion Bf. 



(5.8) 



Ti((a,q),(t,r))= [ (v[B a , s < q - x]P[B a 2 t > r - x] 

— P[B a 2 s < q — x, B a 2 t > r — x]j dx 
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and 

/•oo 

r 2 ((»,«),(*,r)) = / P[-B„i, <5-i]P[B„a,<r-i]<ic 

(5.9) * 

+ / P[-B CT 2 S > <? - x]P[B CT 2 t > r - x] dx. 

J —oo 

By Lemma 5.2, the desired limit (5.1) follows from showing 

(5.10) E w {e iW ™) e~ CT 9 /2 in P-probability as n ^ oo. 

This limit will be achieved by showing that the usual conditions of the Lindeberg-Feller 
theorem hold in P-probability: 

(5.H) y, -+ a ° 

\m\<a(n)y/n 

and 

(5.12) Y, E"(\U rn \ 2 l{\U m \ >e}) -+0. 

\m\<a(n)y/n 

The standard Lindeberg-Feller theorem can then be applied to subsequences. The limits 

(5.11) -(5.12) in P-probability imply that every subsequence has a further subsequence along 
which these limits hold for P- almost every u. Thus along this further subsequence W* 
converges weakly to M(0, a 2 .) under P^ for P-almost every to. So, every subsequence has a 
further subsequence along which the limit (5.10) holds for P-almost every oj. This implies 
the limit (5.10) in P-probability. 

We check the negligibility condition (5.12) in the L 1 sense. 

Lemma 5.3. Under assumption (1.19), 

(5.13) lim V E\\U m \ 2 l{\U m \ > e} 
Proof. First 

/ N V 

Uli = ( ^0i[t/ m (^)l{m > 0} - V m (U,ri)l{m < 0}] J 
^ i=i ' 

N N 
< CY U m{%nfl{m > 0} + CY V m(ti,r l ) 2 l{m < 0}. 
i=i i=i 

The arguments for the terms above are the same. So take a term from the first sum, let 
(t,r) = (tj,rj), and the task is now 

a(n)y/n 

(5.14) lim V E[U m (t,r) 2 l{\U m \ > e}] = 0. 

m=0 

Since 

\U m \ < Cn- l / A [ m {m) +E w {r 1 {m))\ 



n— >oo 

\m\<a(n)^/n 
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and by adjusting e, limit (5.14) follows if we can show the limit for these sums: 

a(n) y/n 

E[U m (t,r) 2 l{r, (m) > n^e}) 

m=0 

(5.15) 

+ £ E[U m (t,r) 2 l{E"(ri {m))>n l l A e}]. 

Abbreviate 

A m = {X™< [nvt\+rJR}. 

The terms of the second sum in (5.15) develop as follows, using (5.5), the independence of 
cj—oo _i and (Do,ooi and the shift invariance: 

E[E"(U m (t,r) 2 )l{E"( m (m)) > n 1 /^}] 

< n- l l 2 E[{E"( m (m)) + Yax"{ m {m)))l{E"{ m {m)) > n l ^£}]P{A m ) 

= n- 1 / 2 E[(^( % (0)) + Var-(r7o(0))) 1(^(770(0)) > n l ' A e}] P(A m ). 
Since the averaged walk is a walk with bounded i.i.d. steps, 

a(n)y/n 

(5.16) Yl P ( Am ) ^ E [( X nt ~ L n ^J " rV^)~] < C{n 1 ' 2 + 1). 

m=0 

Thus 

a(n)y/n 

ElU^rfliE^voim^yn^e}] 

m=0 

< CE[(E»( m (0)) + Var^r/o^l^MO)) > n l ' A e}] . 

The last line vanishes as n — > oo by dominated convergence, by assumption (1.19). 

For the first sum in (5.15) first take quenched expectation of the walks while conditioning 
on rjo, to get the bound 

^o^m(t,r) 2 ]<2n- 1 / 2 P-( J 4 m )[r ? o(m) 2 + ^(r ?0 (m)) 2 ]. 

Using again the independence of oj-oc-i and u)q )OQ , shift-invariance, and (5.16), 

Y E[U m {t,r) 2 l{ m (m)>n l l A e}] 

m=0 

a(n)y/n 

<Cn^ 2 Y ^ m )-^[(r/o(0) 2 + ^(r ? o(0)) 2 )l{r /0 (0)>n 1 / 4 e }] 

m=0 

< ^[(^(O) 2 + E"( m (0)) 2 )l{ Vo (0) > n^e}] 
The last line vanishes as n — > oo by dominated convergence, by assumption (1.19). □ 
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We turn to checking (5.11). 

E E ^m 

\m\<a(n)y / n 

= E E [ 1 {m>o}E J {U m (ti,ri)U m (t j ,r j )) 

l<i,j<N \m\<a(n)yfii 

+ l{m<0}^(K i (^,^)y m (t i ,r i )) 

Each quenched expectation of a product of two mean zero random variables is handled in 
the manner of (5.5) that we demonstrate with the second expectation: 



,r}o(m) 

-1 10. ^, ,,, I 

n 



r) (m) 



l{X% k > [nv ti \ + r iy fr}, E ^ntf > MiJ +fiV^}) 
fc=i £=i ' 

+ n- 1 /2 V ar-(r ? o(m))P-(X- > |mrfij + r^P^X™ > [nvtjj +r jy fc). 
After some rearranging of the resulting probabilities, we arrive at 

E 

\m\<a(n) y/n 



n 



-1/2 



Y, Mi E ^(%H) 



l<i,i<JV |m|<a(n)V^ 

x {P^(X™ < [ntrfij +r jV ^)P-(X- > [nvtjj +r^) 

\m\<a(n)y/n 



(5.17) 



+ l {m < 0} P"(X™ > [nvU\ +r i ^)P^{XZ j > [nvtjj +rjVn-)} 

The terms above have been arranged so that the sums match up with the integrals in 
(5.7)-(5.9). Limit (5.11) is now proved by showing that, term by term, the sums above 
converge to the integrals. In each case the argument is the same. We illustrate with the 
sum of the first term with the factor Var w (r/o(?n)) in front. To simplify notation we let 



30 M. JOSEPH, F. RASSOUL-AGHA, AND T. SEPPALAINEN 

((s, q), (t, r)) = ((ti,ri), (tj,rj)). In other words, we show this convergence in P-probability: 
S (n) = n~ l > 2 Var w (7 ?0 (m))P w (X™ < [nvs\ + q^i) 

0<m<a(n)y/n 

(5.18) x P W PC < [nvt\ +ryfti) 



— > 0", 

n—>oo 



poo 

2 / P[B a 2 s <q-x]P[B a 2 t <r-x}dx = I. 
Jo 



The proof of So(n) — > I is divided into two lemmas. Let 
5i(n) =n~ 1/2 ^ Var aJ (? ?0 (m)) 

/g 0<m<a(n)v^i 



xP(Ba<?-t PLEU<r 



771 \ / _ 771 

^)p(s^ t <r-- p 

77/ V -v/71 



Lemma 5.4. lim E\S (n) - Si(n)\ = 0. 

n— >oo 

Proof. By the quenched central limit theorem for space-time RWRE [11], for each x G R 
the limit 

holds for P-a.e. w. Since these are distribution functions (monotone and between and 1) 
with a continuous limit the convergence is uniform in x. Set 

D n (u) = sup ^{X^ < [nvs\+x^i)P UJ {X nt < [nvt\ + y^/Ti) 

- P(B a 2 s < x)P(B a 2 t < y) | 

and then D n (oj) — > P-a.s. By shift-invariance 

E|5 (n) - Si(n)| < n" 1 / 2 ^ E Var T ° ' " l w ( 770 ( ) ) 

0<m<a(n) y/n 

P To '- mOJ (X ns < [nvs\ + qy/n - m)P T °' mUi (X nt < [nvt\ +ry/n-m) 



= P (B a 2 t <r- — 
n / v V77 



P I -B CT 2 S < a t= P B a 2 t < r 



<rT x l' 1 E[Var T «- w (7 ?0 (0))Z)„(ro, m a;)] 

0<m<a(n)y/n 

(5.20) < 2a(n)E[Var w (? ? o(0))A i (u;)]. 

Moment assumption (1.19) and dominated convergence guarantee that E[Var w (7/o(0))Z? n (a;) 
0. Thus we can take 

(5.21) a(n) = ( sup E [Var w (j7o(0))D fc ( w )] N 



k:k>n 

to have a(n) oo while still line (5.20) vanishes as n — > oo. □ 

The choice of a(n) made above depends on s, t but that is not problematic since we have 
only finitely many time points tj to handle. 
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Lemma 5.5. lim ElS'i(n) — II = 0. 

n— >oo 

Proof. First we discard tails of the sum and integral. Given e > 0, we can choose a large 
enough c < oo such that 

S*(n) =n~ 1 ' 2 Var w (r ?0 (m)) 

0<m<c v / n 

x P(B„ 2j , < a - — ]P(-B CT 2 t < r 



satisfies E|5i(n) — S*(n)\ < e, and so that 



r = o%[ P[B a 2 s <q-x]P[B aH <r-x]dx 
Jo 



satisfies I — I* < e. Thus it suffices to prove S*(n) — > I* . 

Next, since the Gaussian distribution functions are Lipschitz continuous, 



S*(n) - I* = n- 1 / 2 ]T [Vat^C^M) - a 2 ] 

0<m<Cy/n 

* p ( B -^-^) p (^-^)+o(„-^ 



Introduce an intermediate scale 1 << L « y/n and use again Lipschitz continuity of the 
probabilities: 



Sl(n)-r = -^ E (l E Vai"fo,(m)) - A 

0<j<£^-l m=jL+l 

* { p ( s - * « - i m** * * - i ) + °(^) } + § + °<"- i/2 ) 



The error term i? n consists of order L terms bounded by | Var"' (770(777.)) — Oq\ that appear 
because the collection of summation intervals (jL, (j + 1)L] may not exactly cover the 
original summation interval < m < Cy/n. It satisfies Ei2 n < CL. Finally, bounding the 
probabilities crudely by 1 and by shift-invar iance, 



E\S{(n)-I*\ < CE 



1 L 

£Var"Mm))-a 2 



L 



+ 0{Ln~ 1 / 2 ). 



This vanishes as we let first n — > oo and then L — > oo and apply the L 1 ergodic theorem. □ 

Limit (5.18) has now been verified. All terms in (5.17) are treated the same way to show 
that they converge, in L X (P) and therefore in P-probability, to the corresponding integrals 
in (5.7)-(5.9). This verifies limit (5.11). Since both (5.11) and (5.12) have been checked, 
the Gaussian limit in (5.10) has been proved, as explained in the paragraph following (5.12). 
The proof of Proposition 5.1 and thereby also the proof of Theorem 1.4 are complete. □ 
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